Relaxation in time-dependent current-density functional theory 
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We apply the time-dependent current-density functional theory to the study of the relaxation of 
a closed many-electron system evolving from an non-equilibrium initial state. We show that the 
self-consistent unitary time evolution generated by the time dependent exchange-correlation vector 
potential irreversibly drives the system to equilibrium. We also show that the energy dissipated in 
the Kohn-Sham system is related to the entropy production in the real system. 
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The irreversible time evolution of a many-electron sys- 
tem which is initially prepared in a non-equilibrium state 
is the subject of strong research interestllj. The system 
might be, for example, a large molecule J2J or the electron 
liquid in a quantum well, and the initial non-equilibrium 
state might be created by laser excitation, or by the ap- 
plication and subsequent removal of an electric field [2j . 
From experience we know that even if the system is closed 
and isolated, i.e. no exchanges of particles and energy are 
allowed, the interaction between the particles drives the 
system towards the state of maximum entropy, i.e., the 
equilibrium state. When this happens we say that the 
system has relaxed At each stage of the relaxation 
process the system has a certain amount of "mechanical 
energy" which could in principle be channeled into work. 
One of the tasks of theory is to develop efficient tools -in 
alternative to the solution of the many-body Schrodinger 
equation- to calculate the time evolution of the mechan- 
ical energy. One would like to do this starting from first 
principles, yet without knowing the exact state of the 
many-body system as a function of time. 

Various formulations of the problem of relaxation of 
a quantum system have been proposed over the years. 
Most of these formulations couple the system of inter- 
est to some kind of external environment -the so-called 
thermal bath 0. When the information about the mi- 
croscopic state of the thermal bath is discarded this ap- 
proach leads to an effective dynamics of the system that is 
no longer unitary. The main weakness of this approach 
lies in the fact that the thermal bath and its coupling 
to the system must be modeled in an essentially ad-hoc 
manner. In this paper we will describe an alternative 
approach, based on the time-dependent current-density 
functional theory (TDCDFT) [E El 0, IE 0, which al- 
lows a first-principle treatment of relaxation in terms of 
a self-consistent unitary time evolution . 



We assume that the external potential acting on the 
system does not depend on time and that the system is 
closed and isolated. We thus consider the time evolution 
of a many-electron system described by the Hamiltonian 



H 



^ 2m 



W 



dr V(r)n(r), 



(1) 



where h(r) is the density operator, V(r) is the external 
potential and W describes the electron-electron interac- 
tion (here and in the following h = 1). A possible way 
to generate the initial non-equilibrium state \ipi) is by 
applying an external potential at t = — oo and allow- 
ing the system to relax in the presence of this potential. 
The external potential is then switched off at the initial 
time t = and the system evolves henceforth according 
to the Schrodinger equation id t \^{t)) = H\^>(t)). Due 
to the presence of the electron-electron interaction the 
system will relax to an equilibrium state of H. Because 
the system is closed and isolated during its time evolu- 
tion, the energy is conserved and equal to its initial value 
Ei = (tpi\H\ipi). This implies that the final equilibrium 
state will not be the ground state of H, but the most 
probable state with total energy Ei. In this final state 
the entropy, as well as the temperature will have finite 
values. Our task is to calculate the time evolution of the 
mechanical energy between the initial and the final state 
of this evolution. 

The central tenet of TDCDFT is that the correct time 
evolution of the particle and current densities of the 
many-electron system can be obtained from the time 
evolution of a non-interacting reference system, known 
as the Kohn-Sham (KS) system, described by the state 
\^Ks(t)). This state evolves from an initial state \^KS,i) 
-constructed to give the initial density and energy of the 
system- according to the equation 



idt\*Ks(t)) = H KS (t)\^Ks{t)) 

where the KS Hamiltonian is 
H KS (t) = ^^fe + eA cc (f J ,i)] 2 



+ / dr h(r)V hxc (r,t) + / dr n(r)V(r) 



(2) 



(3) 



Vhxc{ r ,t) is an effective scalar potential that includes 
the Hartree potential (h) as well as the ground-state 
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exchange-correlation (xc) potential of static density func- 
tional theory (DFT) evaluated at the instantaneous den- 
sity n(r,t). A xc (r,t) is the xc vector potential, which 
is a functional of the density and of the time-dependent 
current density j (r, t) , and whose role is to enforce the 
correct time evolution of the current. The time depen- 
dence of Hks arises from its being a functional of n and 
j. In practice, both Vh X c( r , t) and A xc (r, t) need to be ap- 
proximated, but these approximations are non-empirical, 
in the sense that they are based on exact constraints com- 
ing from the physics of the homogeneous electron liquid 
y|- For Vhxc(r,t) the simplest and most popular approx- 
imation is the local density approximation (LDA) , which 
expresses the xc potential as a functional of the density 
in the following manner: 



(4) 



Vhxc(r,t) = / dr' -n(r',t) + n xc (n(r,t)) 

J \r-r'\ 



where fJ, xc (n) — de xc {n) j dn and e^in) is the xc energy 
density of a uniform electron gas [j| . 

As for A xc , a quasi-local approximation in terms of 
the current density was derived bv_ Vignale and Kohn 
llOl. and can be cast in the forr 



e-r-^zc k(r, t) = } , ViCjfJr, t) (5) 

dt ' v ' n(r,t) xcK ' ; w 

where i, k denote cartesian components and the stress 
tensor a xc is given by 



af c ir,t) =C(r,t)6 ik V-v 



+ r)(r, t) ( V 4 w fe + Vfcw 4 - V • v 



(6) 



where v{r, t) = j(r, t)/n{r, t) is the velocity, and we have, 
for simplicity, neglected retardation effects in the rela- 
tion between a xc and the velocity. The quantities r\ and 
C are the visco-elastic constants of the homogeneous elec- 
tron liquid and can be calculated from the linear response 
functions of the latter , evaluated at the instantaneous 
density n(r, t). 

We want to point out that the Eqs. define a 

self-consistent unitary evolution which is designed to ap- 
proximate the true density and current of the system. It 
is the dependence of the xc vector potential on the ve- 
locity that breaks the time-reversal invariance of the KS 
hamiltonian and allows for the possibility of irreversible 
relaxation. 

We now define the key quantity of this paper, the 
"Kohn-Sham energy" 



E{t) =(^Ks(t)\H KS (t)\^ KS (t)) 



dr V hxc (r,t)n(r,t) + E hxc [n(t)} 



(7) 



where E} lxc is the Hartree + xc energy functional of the 
ground-state DFT. Eh xc and Vh xc are related by the iden- 
tity Vhxc = SEh xc (n)/Sn. In static DFT, E{t) is just the 



familiar expression from which the ground state energy 
is calculated. In TDCDFT, E(t) is not the true energy 
of the interacting many-particle system, but, as we will 
show below, it decreases monotonically with time going 
from Ei at t = to some final value E / > Eq at t = +oo, 
E being the true ground-state energy of the system. 

To prove this point we observe that a straightforward 
calculation shows that (we used SHks I '$A XC = j) 



dE f i ( \ 
— =Jdrj(r,t) 



dA^ 



dt 



Eq. JSJ immediately allows us to identify dE/dt as the 
work done on the system by the self-consistent xc vector 
potential. By putting Eq. JSJ in Eq. we get 

^ = -2 J dr 7/Tr (u - ^Trc^ - J dr C(Trc) 2 , (9) 



where, uj 



(Vi«j 



VjVi)/2. The non-positivity of 
dE(t)/dt follows from the fact that the viscosity con- 
stants 77 and £ are positive 0- 

Next we prove that E(t) — > Ef > E for t — > 00. We 
begin by observing that, if t Axc = J2i(P* + &A xc ^_j2m is 
the KS kinetic operator, then at a fixed time t |l8| . 

E(t) > min (^(t)\f Axa{t) +V\y(t))+E hxc [n(t)} 

> E [A xc (t), V] + E hxc [n(t)} - E hxc [n{t), A xc (t)\ 

(10) 

where Eq [A, V] is the instantaneous ground state energy 
of the system, in the presence of static potentials V and 
A. Ehxc[n,A] is the Hartree+xc energy functional of the 
system in the presence of a static vector potential A(r). 
Now it is known ^3] that in LDA and up to second-order 
in B(r) = V x A(r), one has E hxc [n] - E hxc [n,A] = 
J dr c(n(r))B 2 (r) > because c(n) (— one half the dif- 
ference between the orbital magnetic susceptibilities of 
the interacting and non-interacting homogeneous elec- 
tron liquid |l2() is a positive quantity. Assuming that 
this inequality is generally true we can then neglect the 
last two terms on the right hand side of Eq. (|10fl to 
arrive at E{t) > Eo[A xc (t),V}. We now observe that, 
because the electron-electron Coulomb interaction is a 
positive operator, the ground-state energy of an inter- 
acting many-particle system is always greater than the 
ground-state energy of the same system where electron- 
electron interaction is turned off. We then conclude 
that E Q [A XC ,V] > E^[A XC ,V] where E^ is the instan- 
taneous ground state energy of the non-interacting sys- 
tem. On the other hand, the ground state energy of 
a non-interacting fermion system is always greater than 
the ground state energy, E^' l b , of a non- interacting boson 
system in the same external potentials. We then have 



E{t) > E™[A XC ,V] > E™[A XC1 V] > E™[A X 



0,V], 

(11) 
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FIG. 1: a) Plot of the dipole moment d(t) — J dzzn(z,t). 
b) Plot of the Kohn-Sham energy E(t) - E . c) Plot of 
AS(t)/ s /2~/(Ei - Eoj as given by Eqs. JTBJ and JTJ. A fit 
with E(t) in Eq. of the data in plot b) gives F = 0.0034 
a.u. We have always used £ — 0.01 mV/nm. 



where we have used in the last step the diamagnetic in- 
equality, which states that the ground state energy of a 
spinlcss boson system in the presence of a vector poten- 
tial is always greater than the ground state energy of the 
same system without the vector potential ^jj. Because 
Eq\[A xc = 0, V] is a time-independent quantity we have 
arrived at our conclusion that E(t) is bounded. Eqs. © 
and (|llfl allow us to conclude that the limit t — > oo of 
E(t) is finite and that lim^oo dE/dt = 0. From Eq. © 
and the requirement that j vanishes at infinity, we de- 
duce that we can have dE /dt = if and only if V • v = 
and V x v — 0, which imply v = everywhere (we con- 
sider only simply connected geometries). In turn this 
implies that the current j is zero everywhere, and there- 
fore, through the continuity equation, the density ap- 
proaches a stationary limit rif{r). This can only happen 
if the asymptotic state is one of the eigenstates of the KS 
Hamiltonian with density nj{r). In the following, we will 
disregard the unlikely cases when the system gets stuck 
in some excited state, and we will assume nj{r) = no(r), 
where no(r) is the ground-state density. The time evo- 
lution of Eq. @ will then drive the KS system to the 
ground-state of the corresponding KS hamiltonian. 

To demonstrate the above ideas in a concrete exam- 
ple, we have computed numerically the time evolution of 
the electronic dipole moment d(t) = J dzzn(z,t) (Fig.^ 
panel a) and the energy E(t) — E (Fig. ^ panel b) in 
a one-dimensional quantum well, such that the electrons 
are confined in the z-direction, and yet are free to move 
in the x — y plane [l4| . The single-particle states for this 
system are of the form ip n (z)e lk ' r where k and r are vec- 
tors in the x — y plane. The system is initially prepared 
in a state that corresponds to the KS ground-state in the 
presence of a uniform electric field £ = 0.01 mV/nm in 
the z-direction. This state is homogeneous in the x — y 
plane, and this property is maintained throughout the 
subsequent time evolution. This means that all the elec- 




FIG. 2: The motion in the z-direction is coupled through 
the Coulomb interaction to a 2DEG in the x — y plane. The 
energy is lost to the motion in the z-direction at a rate given 
by dE/dt. 

trons share the same state of motion in the z-direction 
-a time-dependent state that can be written as a linear 
superposition of the lowest-lying stationary states ip n (z) 
(typically the two lowest ones suffice) . The xc vector po- 
tential A xc has been approximated in the form of Eq. @ 
[l9|. We see from Fig. 1 that the condition dE/dt < 
is verified at all times. Moreover the energy E(t) is well 
approximated by the analytic form 

E(t)-E ^(E t -E )e- rt 1 (12) 

where the parameter T can be estimated from the slope 
of the curve in Fig. 

We now wish to argue that E(t) — Eq is the maximum 
work that can be extracted from the system at a given 
time, and that its time derivative is related to the rate of 
entropy production 0]. To do this, let us take a closer 
look to our quantum well system (see Fig- EJ- In this sys- 
tem we can recognize a macroscopic degree of freedom 
-the z-dependent wave function in which all the elec- 
trons reside- and a "thermal bath" -the two-dimensional 
electron-gas (2DEG) in the x — y plane. The coupling be- 
tween the two is induced by the Coulomb interaction: no 
empirical modeling is needed. As the macroscopic state 
of motion in the z-direction evolves with time, energy is 
transferred from this motion into low-lying excitations of 
the 2DEG in the x—y plane. These excitations are, in the 
simplest case, double electron-hole pairs with zero total 
momentum. Let us study in detail this energy exchange. 

We start from an initial state in which the 2DEG is 
in the ground-state but all the electrons occupy an ex- 
cited state of motion in the z-direction. At the end of the 
relaxation all the electrons are in lowest energy state in 
the z-direction, but the 2DEG is in an excited state. The 
total energy is conserved and given by Ei. On the other 
hand, the infinite-time limit of E(t) is the ground-state 
energy of the system, Eq, as one can see by noting that 
the final Kohn-Sham hamiltonian, according to Eq. (0) 
coincides with the Kohn-Sham hamiltonian of static DFT 
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which, by definition, yields the correct ground-state en- 
ergy. Since the electrons are finally at rest in the z- 
direction, the difference Ei—E must be entirely ascribed 
to the excitation of the 2DEG. This means that Ei — E a 
is the total heat, Q, transferred from the macroscopic 
motion in the z-direction to the 2DEG. 

Denoting by Ti and Tf the initial and final temperature 
of the 2DEG we have 



Q 



dT C V (T) 



(13) 



Ti 



where CV(T) is the heat capacity of the 2DEG. At low 
temperature we have CV(T) = 7T with 7 = TrAm*kB/3 
where A is the area occupied by the gas in the x — y plane 
and m* is the particle mass renormalized by electron- 
electron interactions. Eq. (|13|l provides a quantitative 
relation between the increase in temperature of the "ther- 
mal bath" and the energy that is lost to the macro- 
scopic motion in the z-direction. Because in our model 
calculation Ti = we get for the final temperature 
T f = sJlfiEi - E ) and by recalling that C v /T = dS/dT 
we get the increase in entropy AS = jTf. 

The fact that the energy Ei — Eo is converted into heat 
that increases the temperature and entropy of the 2DEG 
allows us to identify Ei — Eq as the maximum energy that 
could in principle be extracted from the motion in the z- 
direction to do work [2(|. In other words, Ei — E is the 
work we would extract from the system if we could pre- 
vent excitation of the 2DEG during the relaxation of the 
motion in the z-direction. Any excitation of the 2DEG 
will necessarily result in less work obtained from the sys- 
tem. 

We can also use Eq. (|13fl to define a suitable "tem- 
perature" for non-equilibrium states. We begin with the 
observation that Eq. (|13J) can be cast in the form 



dt 



= 0, 



(14) 



which simply states the global conservation of energy. 
This suggest that we define the non-equilibrium temper- 
ature of the system in such a way that the equation 



dE = dS 
dt K ' dt 



is satisfied at all t. The possibility of this definition 
stems from the fact that we can always follow a non- 
equilibrium thermodynamical transformation through a 
series of small quasi-equilibrium transformations. Then 
Eq. lfTK|l is sufficient for the Eq. fHI to be satisfied. 
Eq. (|15|l is in fact a differential equation [2^] for the 
temperature that can be easily solved to give T(t) = 

l[Ei - E(t)] (we have used again C V {T) = -fT for a 



(16) 



2DEG) and finally (see Fig. □ panel c) 
AS(t) = y/2-ylEi - E(t)]. 



This is the desired connection between E{t) and AS{t). 

In conclusion we have shown that the time-dependent 
current-density functional theory naturally opens the 
possibility of a first-principle description of relaxation in 
many-body systems. To show this we have identified a 
quantity, the Kohn-Sham energy, that defines a "time 
arrow" in the evolution of the system. Indeed we have 
proved that this energy monotically decreases with time 
until an equilibrium state is reached. We are also able to 
relate the energy dissipated to the increase in entropy and 
temperature. This thermodynamical approach allows us 
to identify the Kohn-Sham energy with the maximum 
amount of work that can be extracted from the system. 
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